Superintegrability on the two dimensional hyperboloid II 



E. G. Kalnins 

Department of Mathematics and Statistics, 
University of Waikato, 
Hamilton, New Zealand. 
W. Miller, Jr. 
School of Mathematics, University of Minnesota, 
Minneapolis, Minnesota, 55455, U.S.A. 
Ye. M. Hakobyan and G. S. Pogosyan 
Laboratory of Theoretical Physics, 
Joint Institute for Nuclear Research, 
Dubna, Moscow Region 141980, Russia 

February 1, 2008 

Abstract 

This work is devoted to the investigation of the quantum mechanical systems on the two 
dimensional hyperboloid which admit separation of variables in at least two coordinate sys- 
tems. Here we consider two potentials introduced in a paper of C.P.Boyer, E.G. Kalnins and 
P.Winternitz, which haven't yet been studied. We give an example of an interbasis expansion 
and work out the structure of the quadratic algebra generated by the integrals of motion. 

1 Introduction 

Superintegrable systems on the two-dimensional hyperboloid were introduced and developed 
in the papers |T|, 0, |[]. In distinction to the cases of two-dimensional Euclidean space and 
the two-sphere, the classification of superintegrable systems on the hyperboloid is difficult. 
To date only the four potentials studied in || and two more listed in || are known. In 
the present paper two potentials are considered, which were constructed in the work |IJ] but 
have not previously been investigated. These potentials both have only a finite number of 
bound states. At this point we have treated all the potentials that arise by restriction from 
hermitean hyperbolic space. We follow the approach of 0, which contains an introduction 
and motivation. 

The two dimensional hyperboloid is characterized via the cartesian coordinates uq, u\, u>2 
where uj\ — uj\ — uj\ = 1,uj > 1. The requirement u> > 1 means that we consider only 



upper sheet of the double-sheet hyperboloid. Throughout this paper we will consider the 
Schrodinger equation on the hyperboloid in the form (h = m = 1) 

m= (~A LB + v)y = EV (1) 

where V is a potential function and the Laplace-Beltrami operator Alb is written as 

A LB = Kl + K 2 2 -Ml (2) 
Here K 3 , K 2 , Mi generate the Lie algebra so(2, 1) [|, || 

K 3 = uj d Ul + uJid U0 , K 2 = wo^wa + ^2^, Mi = uid U2 - u 2 d ui . (3) 

and 

[K 3 ,K 2 ]=M 1 , [K 2 ,M 1 ] = -K 3 , [K 3 ,M 1 ] = K 2 (4) 

The Schrodinger equation ([!]) for V = separates in nine coordinate systems ||. Introduc- 
tion of a potential breaks the symmetry and, in general, reduces the number of coordinate 
systems permitting separability, usually to zero. We consider the following two potentials 
(see Table), constructed in [Q, for which ([I]) is superintegrable. 

Table 



Potential V(u) 


Coordinate system 


y — a 2 7 2 1 g2 t^o+^l 


Equidistant 
Elliptic-parabolic 

Hyperbolic-parabolic 
Horicyclic 


\r — a 2 1 -,2 uQui 

^-^ + 7 

1 / 2 «2N wg-w? 


Equidistant 
Semi- Hyperbolic 



Recall that ([]]) is superintegrable for a given potential V if it is separable simultane- 
ously in at least two coordinate systems. 



2 First Potential 

The first considered potential is 



a 2 7 2 n2 u + loi 



Vr = - 2 - ( ' \, (5) 
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where a, (3, 7 are positive constants . The corresponding Schrodinger equation admits separa- 
ble solutions in four coordinate systems: equidistant, elliptic-parabolic, hyperbolic-parabolic 
and horicyclic. 



2.1 Solutions of the Schrodinger equation 



1.1 Equidistant coordinates. In this coordinate system 

u>o = cosh 7~i cosh T2, 0J\ — cosh 7~i sinh T2, cl> 2 = sinh T\ 
[ti, r 2 G (—00, 00)] the potential V\ has the form 



^(n,^) 



a 2 1 /3 2 — 7 2 (coshr 2 — sinh r 2 ) 2 



sinh ri cosh T\ (cosh r 2 — sinh r 2 ) 4 



After putting 



*(n,r 2 ) = (coshn)- 1 / 2 ^^)^^) 

we come to the system of equations: 



<PS 2 
dr! 

d 2 S l 



+ 



+ 



-fi 2 - 2/3 2 e 4r2 + 2 7 2 e 2r2 



(2^ " 7) + 



5 2 = 
2a 2 



4 7 cosh 2 7~i sinh 2 ri 



Si = 



(6) 
(7) 

(8) 
(9) 



where \i is the equidistant separation constant. The first equation (^) could be considered as 
a one dimensional Shrodinger equation for the Morse potential || and the orthonormalised 
solution is given by the expression: 



S 2 (t 2 )=S^\z) 



2^r(m + fM + 1) , /2 
^ m!r 2 (// + l) 6 * 2 l( m '^ + i ' 



2/im! 



_ e --/V /2 L^), z = ^2/3e 



2r 2 



V T(m + /i + 1 

where -^mW are the Laguerre polynomials @. The separation constant is quantized as 



(10) 



-2m - 1 



7' 



>/2f}' 



< m < 



2 Vv/2/3 , 



(11) 



The second equation (||) represents the modified Poschl-Teller equation || ||. The orthonor- 
malised wave function is given by: 



S 1 (t 1 ) = S^in) 



2(ji - J2a 2 + 1/4 -In- l)r(/x - n)n\ 



\ T(/i - pa 2 + 1/4 - n)r(l + n + v /2a 2 + 1/4) 
(sinhr 1 )Hv /2 « 2 +V4 (coshri) i-M P y 2 « 2 + 1 /4,-^) (cosh2ri)) 



(12) 



with n = 0, 1, . . . 



quantized energy is 



^ (u — 1 — \j2a 2 + |) , where P^ a '^(x) is the Jacobi polynomial 



E 



N 



--(//- pa 2 + l/A -2n-l) 2 



1 



where iV = m + n is the principal quantum number and the bound states occur for 



< N < 



[2 \V2f3 



2a 2 + 1/4-2 



The orthonormalized total wave function *& nm (Ti,T2) is given by (|7|), fllOD and (|T^) 
The symmetry operator describing this coordinate system is 



Kl - 2(3 



2 f ^O + uh 

\LU — LUi 



+ 2 7 



2 UJ + LOi 



UJ - U) X 



l 2 V 

-2m - 1 + j ^ 

rim 



1.2 Horicyclic coordinates. In the horicyclic coordinates 



x 2 + y 2 + 1 x +y — 1 x 

^0 = ~ , ^1 = ^ , ^2 = — , 



[y > 0, x G (— oo, oo)] the potential V\ is 



a 



+ (3 2 (x 2 + y 2 )-Y 



and the Schrodinger equation has the following form 
1 



~2 V 



Via putting 



d 2 2n 2 d 2 

h-^r- 2(32x2 + h~ 2 ^ 2 + 2 ^ 2 

aar x A oy z 



V(x,y) = i)i{x)i) 2 {y) 



V(x,y) = EV(x,y) 



it admits a separation 



dx 

d 2 j, 2 
dy 2 



i + 2 



+ 2 



7 2 (A 1 + 1)-/3V-^ 
ar 

E 



7 2 (A 2 -1)-/5V + 



^i=0 

^2 = 



where Ai and A2 are the horicyclic separation constants with the relation: Ai + A2 = 1. 
The orthonormalized solutions of the equations (20), fl2"T|) for (-2E + 1/4) > are 



\ T(n 1 + .\/2a 2 + 1/4 + 1) 



x 



e"^(VV2/?a; 2 )^ + v2^+i 



28x z 



(22) 



ife(v) = ^S ,/3) (?/) 



n 2 !( v / 2/3) 1 /2 



\ r(n 2 + ^/-2£ + 1/4 + 1) 



x e-^( V / ^)^v^^LV r ^(^ 
The separation constants A 1; A 2 are quantized as: 



y 2 )- 



Ai 



V2{3, 



V 



(2m + J2a 2 + l/4 + l)-l; A 2 



V 



(23) 



(2n 2 + J-2E + 1/4 + 1) + 1 (24) 



and according to the relation Ai + A 2 = 1, we come to the energy spectrum as in (13). The 
operator characterizing the separation in horicyclic coordinates is: 



{K 2 - M 1 



28 2 u 2 2a 2 (u -u; 1 )> o 2 

h 27 



(uj - Wi) 2 



2^2/3(2^ + ^a 2 + 1/4 + 1) + 2 7 2 



*ni« 2 (s,!/) 



(25) 



i.5 Elliptic-parabolic coordinates. In this coordinate system 

sinh 2 a — sin 2 9 



cosh 2 a + cos 2 6 1 
2 cosh a cos # 



2 cosh a cos 9 



, c<j 2 = tanhatan6>, 



[a > 0, 6* G (— f , |)] the potential Vi has the form: 



Vx(a,9) 



cosh 2 a cos 2 # 
cosh 2 a — cos 2 9 



/? 2 (cosh 2 a sinh 2 a + cos 2 sin 2 



7 2 (cosh 2 a — cos 2 9) + a 2 f 5 1 27: 

Vsinh a sin 6* 



The Shrodinger equation is: 

1 cosh 2 a cos 2 9 



+ 



2 cosh 2 a — cos 2 
d 2 



d 2 2a 2 
23 2 cosh 2 a sinh 2 a + 2 7 2 cosh 2 a 



9a 



sinh a 



2 n 



<9fl 2 



-2^ cos^sin 2 #-2 7 " cos 1 9 



2a 
sin 2 9 



V(a,6) = EV(a,i 



(26) 



(27) 



(28) 
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Putting for the wave function $(a, 6) = S(a)S(9), after separation of variables we get two 
identical equations: 



d 2 5(p) 
dp 2 



+ 



A — 20 1 cosh 2 p sinh 2 p + 2 7 2 cosh 2 p — 



2a 2 



2E 



sinh p cosh p 



S(p) = (29) 



where A is the elliptic-parabolic separation constant and p = a, i6. After changing the 
variables x = cosh 2 p in eq. (|29"D , we obtain 



4x(x-l)^ + 2(2x-l)^ + 
dor dx 



A - 2 / 3 2 x(x - 1) + 2 7 2 x 



2a 2 2E 



x — 1 x 



5 = (30) 



Thus the region x G [1, oo] in eq. ( p0|) belongs to the wave function 5(a) and a; G [0, 1] to 
the wave function S{6). Putting 



Six) = {x- l)Ve- /to /^G(x), 



where 



1 1 „ 1 1 1 / „ 1 

= - + — 1/a 2 + 77> ^ = 7 + — \i-E + 



we get 



d 2 G 1 
Ite 2 + 2 



4 



1 + 4t 1 + 4s 4/5 



4 ^2 



(31) 



(32) 



dx 



x - 1 v/2. 

1 f [2 7 2 - 4/3(1 + 2(t + s))/V2}x + u + V2(3(l + At) + 4(t + s) 2 \ Q _ Q 
4 1 x(a; - 1) J 

If we now substitute 

v 

G(z)=n(z-0i) 



(33) 



(34) 



i=l 



and take into account fl3~2|) , we find that #j satisfies the equation: 



20^1-0, 



^ ^ 1 /3 ^ 
k=i ®k — 9% \/2 



+ 2{i-e i )N + ^-e i 



+ 



7 



'2« 2 + - + l 



0. 



(35) 



The quantization for the energy is given via: 



11 Y 
-2E+- + \ 2a 2 + - + 2N + 2- -!=- 

4 V 4 y/2/3 



(36) 
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and we obtain the expression ([Uf). The separation constant A is: 

v^fe * \V2/3 J ' s/2 
Thus the total solution \I/(a, 9) is represented as: 



(37) 



= S Np (a)S Nq (9) = (sinh a sin 0) W 2a2 +i (cosh a cos ^-^^^^ 
• exp | --^=(cosh 2 a + cos 2 0) j JJ(cosh 2 a- t ) (cos 2 0-0*) (38) 

where p and q is the number of zeroes for the wave functions S(a) and S(9) in the regions 
[0, 1], [1, oo] correspondingly; and the total number of zeroes is N = p + q. 

Eliminating the energy E from equation (p0[), we see that the additional integral of 
motion here is 



1 f d 2 d 2 

— — tt, tq—\ cosh 2 a— -7T + cos 2 9— — - 2 / 3 2 (cosh 4 a sinh 2 a 

cos 2 9 - cosh a I aa 2 o9 2 

+ cos 4 sin 2 0) + 27 2 (cosh 4 a-cos 4 ^) - 2a 2 (coth 2 a - cot 2 0)Wjv M (a, 9) 



(K 2 - M 1 ) 2 - K\ + 20- 



(Wq + Wi) 2 + w\ 
(w - Wt) 2 



2a 2 



w 2 



- 4 7 2 — ^ — \y Npg {a, 9) = X^Npqia, . 



(39) 



i.^ Hyperbolic-parabolic coordinates. In this coordinate system 



cosh 2 6 + cos 2 9 
2 sinh 6 sin 9 



sinh 6 — sin 9 
2 sinh 6 sin 9 



u>2 = coth b cot 0, 



[6 > 0, 9 e (— f , |)] the potential Vi has the form: 



^i(6,0) 



sinh 2 6 sin 2 9 
sinh 2 6 + sin 2 9 

.2/ • u2 7 . • 2 



/5 2 (sinh 2 6 cosh 2 6 + sin 2 cos 2 

2f 1 1 



- 7 (sinh 6 + sin 9) + a ( — — — — t— 

Vcos 2 cosh 2 b 



The Shrodinger equation is 



1 sinh b sin 9 



2 sinh b + sin 9 







2a 2 



— - 2/T sinh 2 b cosh 2 6 + 27 2 sinh 2 b H — 

ob 2 cosh 6 



+ |^-2/3 2 sin 2 0cos 2 + 2 7 2 sin 2 0-^ 
o9 2 cos 2 6* 



#(6,0) = 



(40) 



(41) 



(42) 
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Putting for the wave function ty(b,9) = S(b)S(9), after separation of variables we get two 
identical equations: 



d 2 S(p) 
dp 2 



^ — (3 2 sinh 2 p cosh 2 p + 7 2 sinh 2 p + 



E 



cosh 2 p sinh 2 p 



S{p) = (43) 



where r is the hyperbolic-parabolic separation constant and p = b, id. After changing the 
variables x = sinh 2 p in eq. (fE|), we come to the equation 

2a 2 2E 



4x(x + l)^f + 2(2x + l)^ + 



t - 2(3 2 x(x + 1) + 2Yx + + 

X + 1 X 



Choosing 



N 

P{x) = (1 + ijVe-^^nt^ - Oi) 

i=i 



5 = 0. (44) 



(45) 



where t and s are given by the formulas (0), we obtain the energy spectrum fl3"6"|). Here 9{ 
satisfies the equations 



20,(1 + 0; 



E 



i 



2(1 + 9 i )N+^-Le i + 



r 



y/2/3 



'2a 2 + - - 1 



0. 



The separation constant r is: 



T 



(46) 



(47) 



so the total solution 0) is represented as: 



q Nlk (b,9) = S Nl (b)S Nk (6) = (cosh b cos 0)3+v^+i(siiih6 sin 0)^ v^^+i 



/5 , • 



JV 



exp j--^=(sinh 2 6 - sin 2 0) J JKsinh 2 b - 4 )(sin 2 + f ). 



(48) 



The total number of zeros is iV, and k of them are located in the interval [—1,0] and I are 
in [0, 00]. 

Each solution ^/Nik{b,9) satisfies the eigenvalue equation 



L^ Nlk (b,9) 



I f d 2 d z 

— r- 277 s sinh 2 b-~-r — sin 2 0^^- — 2/3 2 (cosh 2 6 sinh 4 6 



<9 2 



sin b + sin I 



db 2 



dd 2 



- cos 2 sin 4 0) + 2 7 2 (sinh 4 b - sin 4 0) + 2a 2 (tanh 2 b + tan 2 0) J ^ Nlk (b, 1 



(Jf 2 - M^ 2 - Kl + 2 ^K + ^) 2 -^ 2 _ 2ft2 ^0 



- 4 7 2 



U7 - Wl 



^mk(b,9) =TV Nlk (b,9). 



W-2 



(49) 
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2.2 Algebra 

Among the operators {Li, L 2 , L 3 , L 4 }, corresponding to the four separable coordinate sys- 
tems, only two are independent, as 



L3 — —L 2 — L\, L4 — L 2 — L\. 
Consider the operators N±, N 2 and R where 
N x = L u N 2 = L 2 -2 1 2 , 

R = [N 1 ,N 2 ] = 2{K 3 ,{K 2 ,M 1 }}-2{K 3 ,K 2 }-2{K 3 ,M 2 } + 8 

1Q(3 2 lo 2 



(50) 



a 



2 (UJq — UJ\ 
UJ 2 



+ 



+ p 



- 4 



UJ 2 

uj - ui 



2i 



K 3 + 



7 2 + 2a' 



(uj - LOx) 2 
2 



(UJ K 2 - LO X Mx) + 



u - ^1 



(Mi - K 2 ) 



.2 . 0.2/ ~ i"> - ^'i \ * _ 1 + 2 ^2 



UJ 2 



(u Q - UJx) 2 



We have 
[R,N 2 ] 
[R,Nx] 

R 2 



8N 2 + U(3 2 H + 16 7 2 iV 2 + 32/3 2 iVi + 16/3 2 (1 - 4a 2 ) 
-8{iVi, iV 2 } - 32 7 2 # + 16iV 2 - 16 7 2 iVi + 167 2 (1 - 2a 2 ) 

-{iV 2 , iV 2 , iV 2 2 + 32(3 2 N 2 + 128(3 2 H 2 + U^ 2 HN 2 + 128/3 2 #iVi 

3 3 



(51) 

(52) 
(53) 



+ IfrfiNu N 2 } + (™ + 256a 2 ) (3 2 H + (e4a 2 7 2 - ^fl 2 ) N 2 



+ 



352 



128a 2 ) P 2 Nx + (128a 4 /3 2 + 128 7 4 a 2 - — a 2 /3 2 - -/3 2 - 48 7 4 ) 
/ 3 3 



where {A, B} = AB + BA and 

{A, B, C} = ABC + ACB + BCA + BAC + CAB + CBA. 
The integrals of motion Ni, N 2 and H generate a quadratic algebra. 



2.3 Interbasis expansion 

For a fixed value of energy, we can write the equidistant wave function (|7|) in terms of the 
horicyclic ones fll9D as 



ni+n 2 



(54) 



m=0 



where n\ + n 2 = n + m. The connection between the equidistant (a, b) and horicyclic (x, y) 
coordinates is 

x = e fe tanha, y = e b — - — . (55) 

cosh a 



Going over to the horicyclic coordinates in the left side of expansion fl54|) , then considering 
the limit b — > oo and using the asymptotic formula for Laguerre polynomials M 



lim L"(x) 

x—>oo n x ' 



n! 



(56) 



we see that dependence on 6 cancels on both sides of fl54|). Now using the orthogonality con- 
dition for the angular wave functions (|T^) we find the following expression for the interbasis 



coefficients W™™ 2 : 



W nm 



■IT, 



mm 



\y/2P(ji -d-2n- l)T(ji + m + l)r(// 



n] 



\ n x \n 2 \^Y(ni + d+ l)T(n 2 + d+ l)T(n + d + -d-n 



-Bl 



where 



r>nm 

nin 2 



+00 



(sinh a 



,l+2d+2ni 



(cosh a 



, 1— 2/i— 2m 



PV'-A (cosh 2a) da 



(57) 



(58) 



and d = y 2a 2 + 1/4. The integral B™™ n2 can be evaluated by expressing the Jacobi polyno- 
mial through the hypergeometric function 2 F X 0: 



, r(n + /? + !) 
y 1 T((3 + l)n\ 2 1 



-n, n + a + (3 + 1 
+ 1 



1 + x 



(59) 



Representing the function 2 -Pi as a series we come to a sum of integrals, each of which can 
be calculated by using the formula [^]: 



+00 



(sinhr) a (coshr) -/3 dr = -B 



1 + a (3 — a 



[Rea > -1, Re{a - (3) < 0]. (60) 



We thus obtain 



If™ 

mn2 



.iy» 



m!\/2/3(A* - d - 2n - + m)r(n x + d + 1) 
2 \| n!ni!ri 2 !/ir(n 2 + d + l)T(n + d + l)r(/i - n - d) 



(61) 



r(/i)r(/i + m — d — ni — 1) 



-3^2 



-n, n + d — fi + 1, 1 — ix — m 
1 — 11, 2 + ni + d — (i — m 



T(/i — n)T(n + m) 

Alternatively, by using the formula [H| for the Hahn polynomials h^'^(x, N). 



h^\x,N) 



-l) n T{N)T{f3 + n + 1) 
n!r(iV-n)r(/3+ 1) 



■3^2 



— n ; a + (3 + ra + 1 
(3 + 1 ■ 1-N 



x 



(62) 



we obtain the following expression for the expansion coefficients 

(-1)' 



n\U2 



m 



\n\y/2P(n - d - 2n - l)(/i + 



m) 



2 \ ni\n 2 \fiT(n + d + l)T(/i — n — d) 



V[n x + d + l)V{ji-n) 



m) 



\ T(n 2 + d + l)T(fi + 

■ h&'riQi + m + 1,/i + m - d- m - 1) 
in terms of Hahn polynomials. 



r(/i + m — d — ni — n — 1) 



(63) 
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3 Second Potential. 

The second considered potential is 



^2 K + UJ() 2 



H + ^ 2 ) 2 



(64) 



where a, (3 and 7 are positive constants. The corresponding Schrodinger equation admits 
separable solutions in two coordinate systems: equidistant and semi-hyperbolic. 

3.1 Solutions of the Schrodinger equation 



2. 1 Equidistant coordinates. In this coordinate system 

ojq = cosh T\ cosh r 2 , UJ\ — cosh T\ sinh r 2 , uj 2 = sinh T\ 
[ti, t 2 G (—00, 00)], the potential V 2 has the form 

V^ri, r 2 ) = 

After putting 



a 2 I a 2 — /3 2 + 7 2 cosh r 2 sinh r 2 



sinh 2 n cosh 2 n (cosh 2 r 2 + sinh 2 r 2 ) 2 
^(r 1 ,r 2 ) = (coshr 1 )- 1 / 2 5(r 1 )Z(r 2 ) 



we arrive at two equations: 
d 2 S 



dr 2 



dr 2 



+ 



+ 



2£- 



2(a 2 -/5 2 ) +7 2 sinh(2r 2 ) 



1 , /^ 2 -i 



cosh 2 (2r 2 ) 

2a 2 



5 = 



4 cosh 2 n sinh 2 t\ 



S = 



(65) 

(66) 
(67) 

(68) 
(69) 



where \i is the equidistant separation constant. 

Let us consider the first equation ( |68|) . The substituon x = sinh2r 2 transforms this 
equation to 



d 2 S 



dS 



4(1 + x z )— +4x— + 



dx 



dx 



2 2(/j 2 -a 2 )- 7 V 
_/i + (l + a; 2 ) 



5 = 



(70) 



where the physical region is a; G (—00, 00). The equation (|68| ) has three regular singularities 
in the point x = —i, i, 00 and may be solved in term of hypergeometric functions. The 
solution of the equation (|Blf) for a large \x\ can be written as: 



S{X)=A 1 {X-1) 5 4( X + l ) 2 + 4 2 F 1 f ^ tl-l+fi-—— 



+ A 2 (x - i) 



b-ii 1 



x + i)a + 4 2 Fi 



a + 6 + 1 — // 6 — a + 1 — /U 



1 - /i; 



2/ 



2 — X . 



(71) 
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□ 

□ 



□ 



Figure 1: Domain of Convergence 
2, and converges on the circle C with the condition Re(b) < 0. The 

12 



function S(x) exists everywhere inside C except the interval x G (see [I]), since the 

hypergeometric function in ( [73] ) has a cut along the argument G [l,oo)). That means 
that the solution (|73|) along the real axes inside C in general is not a continuous function 
and may have a jump at the point x = 0. Let us now consider the analytic continuation of 
(|73|) inside the circle C 



S(x) 



Al (x — i)2 + 4(x + i) b 2 + 4 



r( M + i)r(-a) 



a+b+l+n 



r( ^a±l ±g)r( -6-a+l+ M )(2i) ^ 



_, / a+b+1 + 11 a + b + 1 - a 

2*1 - . = ; a + 1; 



2 



+ — i) 



2 

l (x + i) 



-+- 

2 ' 4 - 



2i 



-a + b + 1 + fi —a + 6 + 1 — ji 



)(2z) 2 



2 2 
^From equation ( fT2[) follow two possibilities 

a = b* , a = —6* 



-a + 1; 



x 



2i 



(74) 



(75) 



Putting the a = b* (Re(a) = Re(b) < 0) we find that the first term in (f74"|) represents an 
analytic function, while the second term is discontinuous at x — 0. [Note since the both 
terms in equation ( fTif ) transform to each other with replacement a — > —a the choice a = —6* 
means that the first term in ( j74|) is discontinuous while the second term is continous at 
x = 0.] Thus the sufficient condition for the existence of the continuous solution requires 
the relation 

a + a* + 1 " 



so from 



/x + a + a* + 1 
we have 



-2m, 



m = 0,1,2. 



(76) 



H = --zm - i + — \/ z/3 2 - 2a 2 + 1 + ^J(2/3 2 - 2a 2 + l) 2 + 7 4 . (77) 
Finally, the orthonormalized eigenfunction of equation (|68|) may be written in the form 



S(r 2 ) 



(-i)-r(-a), 



'—2m — a — a* — l)T(—m — a* 



\ 7rm!2 a+a * +1 r(— m — a)T(— m — a — a* 
[1 + zsinh2r 2 )^ + 3(l — isinh2r 2 ) 



a I 1 

2 '4 



2 F 1 —m, m + a + a* + 1; a + 1; 



1 + i sinh 2r 2 



\ 



-2m — a — a* — l)m\T(—m — a)T(—m 



7r2 a+a * +1 r(-m - a - a*) 
[1 + i sinh 2r 2 )t+3(l -% sinh 2r 2 )^ + 5P^ a *)(-i sinh2r 2 ) 



(7? 
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where 



2§ 



(2(3 2 - 2a 2 + l) 2 + 7 4 + 2/3 2 - 2a 2 + 1 



+t y v /(2/3 2 - 2a 2 + l) 2 + 7 4 - (2(3 2 - 2a 2 + 1) j . 
The second equation (|69|) is quite like @ and has a solution: 



z( Tl ) = Si^in) 



2{ji - J2a 2 + 1/4 - 2n - l)r(// - n)n! 



x 



^ r(/i - ^a 2 + 1/4 - n)r(l + n + ^a 2 + 1/4) 
(sinh n ) W 2q2 +V* ( cogh ri ( cosh 2ri ) 



with n = 0, 1, .... 

The quantized energy is 



E 



- v/2« 2 + 1/4 - 2n - l) 2 + ~ = -i| 2JV + 2 + y^a 2 + 1/4 



v/2 



2/3 2 - 2a 2 + 1 + v/(2/3 2 - 2a 2 + l) 2 + 7 4 



1 



where N = n + m is the principal quantum number and the bound state occurs for 

o <^ .V < j 



2/? 2 - 2a 2 + 1 + ^/(2/3 2 - 2a 2 + l) 2 + 7 4 - -y/2a 2 + 1/4-1 



The additional operator describing this coordinate system is 



Ki-2(a 2 -P 



2 2 \ ^ 



-2Y 



H + co 2 ,) 2 



1 



2m + 1 - -j= \\2(5 2 - 2a 2 + 1 + y/(2/P - 2a 2 + l) 2 + 7 4 } * nm (r ls r 2 ). 



Semi-hyperbolic coordinates. Here 



(//- e 3 )(i/- e 3 ) 1 _ 1 



2[(e 3 - a) 2 + 6 2 ] 2 26 

(jj, - e 3 )(v - eg) _ 1 _ 1_ 
2[(e 3 -a) 2 + 6 2 ] 2 26 

Q*-e 3 )(i/-e 3 ) 
(e 3 - a) 2 + 6 2 



[( / u-a) 2 + 6 2 ][(z/-a) 2 + 6 2 ]l 1 /2 



(e 3 - a) 2 + 6 2 
[( / u-a) 2 + 6 2 ][(z/-a) 2 + 6 2 ]l 1 /2 



(e 3 - a) 2 + 6 2 
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[v < e 3 < fi], where sinh2/ = (e 3 — a)/b and 2f is the distance between the focii of the 
semi-hyperbolas and the bases of their equidistants 0. 
If we change variables according to 



u 



[Sl + s 2 



-i(s x - s 2 



M) 



the Schrodinger equation becomes 



_d_ 

ds-2 



s —Y + ( s A 

S2 ds 1 ) \ds 3 



S3 



_d_ 

dsi 



S3 



_d_ 



S2 



_d_ 



-E-U* 



2 _ I 

4 



+ 



P _ 1 p _ 1 



<? 2 

■ s 2 



<i 2 
s 3 



^ = 



55) 



with 
Noting 



1. 

4' 



77 



1, 
4' 



2 a 2 ^ 



r7 



^0 ~ ^1 ~ ^2 



2 1 2 1 2 

S l ~r S 2 ~r S 3 



-04 - 7! 

2 V 3 4' 



and considering eq.(|35D we see that the problem we wish to solve using the real coordinates 
ujo,uji and 0J2 is a real case of the corresponding problem on the sphere with coordinates 
Si, S2, S3 and energy e = —E. 

Inverting the relations (|84"D we have 



si 



S2 



S3 = tu 2 - 



Now choose elliptic coordinates on the complex sphere according to 



((j, - ei)(z/ - ei) 
(ei - e 2 )(ei - e 3 ) 



(jm - e 2 ){v - e 2 ) 
(e 2 - ei)(e 2 - e 3 ) ' 



(jm - e 3 )(u - e 3 ) 
(e 3 - e 2 )(e 3 - d)' 



This choice of real coordinates fi, v will work for the real coordinates uik, k = 0, 1, 2 if we 
take ei = e\ = a + a 5 & rea l an d ^ < e 3 < /x. 

In terms of the coordinates /z and ^ the Schrodinger equation has the form: 



(ji-v) 



(/i - e* 2 )(ii - e 2 )(ji - e 3 ) 



-(z/-e;)(z/-e 2 )(z/-e 3 ) 



1 

+ - 



1 1 1 
+ + 



\9* 

dy? ' 2 V yU — e 2 ' /i — e 2 ' /i — e 3 / <9/i 

«9 2 ^ 1/ 1 1 1 \<9# 

<9z/ 2 2\v — e\ v ~ e 2 v — e3/dv 



+ 



ffc2 _ l, (e^ -e 2 )(e* -e 3 ) 2 _ 1 (e 2 - e*)(e 2 - e 3 ) 
M 1 ' ' - 2 4 J (^-e 2 )(^-e 2 ) 



A J (ji-e®(u-eZ) 



2 l A (e 3 - e 2 )(e 3 - 



4 (p-e 3 )(y-e s 
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The separation equations are: 



(p-e 2 )(p-e 2 )(p-e 3 ) 



d 2 ^ 1/1 



+ 



dp 2 2\p — e 2 p — e 2 p — e 3 Jdp 



(P - e^) 



Me 2 -e*)(e 2 -e 3 ) 
4 J (p-e 2 ) 



l, (e 3 -e 2 )(e 3 - e*) 
4 J (p-e 3 ) 



Hp) = 



(87) 



where p — p, v. The operator L 2 with eigenvalue A is 



-4 



(p-u) 



v{p - ei)(p - e 2 )(p - e 3 ) 



<9 2 ^ 1 



1 1 1 \9f 

+ + 



dp? 2 V p — e 2 p — e 2 p — e 3 J dp 



-p[{v - ei)(v - e 2 )(v - e 3 ) 



<9 2 ^ 1/1 



+ t 



+ 



+ 



1 \<9^ 



4 (A*-e5)(f-e5) 



<9z/ 2 2\v — e 2 v — e 2 v — e 3 ) dv . 

{p + v - e 2 ) 



4 (// - e 2 )(z/ - e 2 ) 

(p + is- e 3 ) 



4 - e 3 )(i/ - e 3 ) 



*. (88) 



In order to find the bound state solutions of this system in semi-hyperbolic systems we first 
observe the identity 



+ 



+ 



Qj - e* 2 9j - e 2 6j - e 3 



(uj 2 , — ujf)(6j — a) — 2uj uJib 
(6 j -a) 2 + b 2 

{n-6j){y-6j) 
(9 j -e* 2 )(e j -e 2 )(9 j -e 3 y 



LUn 



Oj - e 3 



(89) 



If we then look for solutions of the form 



3 , 1 N 
* II*/' "MI 



'i 



+ 



s 2 



+ 



1=1 j=l ^Oj - e 2 9 3 ~ e 2 Oj - e 3 

we see that the corresponding zeros satisfy the equations 



(90) 



fci + 1 k 2 + 1 k 3 + 1 



N 

E 



9 m — e* 2 9 m — e 2 9 m — e 3 
For the energy E we have 



Oj) 



0. 



(91) 



E = -\(2N + 2 + k ± + k 2 + k 3 ) 2 + I, 
2 o 



(92) 
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which coincides with the formula (|73|), note (|86|). For the separation constant A we obtain: 
A = -2[k 1 (e 2 + e 3 ) + k 2 (e* 2 + e 3 ) + k 3 (e 2 + e* 2 )] - 2[e 3 k x k 2 + e 2 k x k 3 + e* 2 k 2 k 3 ] 

- 1(4 + e 2 + e 3 ) - 4e 2 e 3 (. 1 + 1) ± - e^ { k 2 + 1) ± 

- 4e 2 e;(fc 3 + l)E 1 (93) 

m=l l^m - e 3 J 

In terms of variables Wi the total wave function is written: 

uj + iui\ kl+ ^ fuj — iui\ k2+ h ^ ^ k3Jr i^j \{uJq — uj\){9j — a) — 2ujQUJ\b u\ 



* = PtPJ F7PJ M* + *n 



a) 2 + 6 2 0, - e 3 



The algebra of second order symmetries for this potential is generated by the operators 
L jk = ( Sj d Sk - s k d S] f + (i - fc|)§ + (i - kl)i (94) 
for j, k — 1, 2, 3 and j 7^ fc. The Hamiltonian of the system is expressed in terms of Ljk as: 

# = + L 13 + L 23 ) - \ £ A; 2 + 2 (95) 

The relevant generators in the real case we are considering are then 

1 2 \ (ojq — iui \ 2 / 1 2 \ /cj + %u)\ x 2 



L 12 = -K'i+ [--k p— ^ + U-* ( 96 ) 



-{M l -iK 2 ) +10 -a --7 - — +a s 

2 V 2 / (uq + ^i) 2 u> 2 



^3 = ^(M 1+ ^ 2 ) 2 +(>- a 2 + V) . A . 2 +a 2 ( "°~r )2 - 08) 
2 V 2 / (lu — luJi) w 2 

The commutation relations and resulting quadratic algebra can then be deduced from the 
relations for the complex forms in terms of the L^. It is easy to show that the additional 
integrals of motion, corresponding to the separation in equidistant and semi-hyperbolic 
coordinates can be written as 

In = -L n +{3 2 -a 2 (99) 

and 

L 2 = e 3 L 12 + e 2 L 13 + e x L 32 - k\(e 2 + e 3 - e x ) - k\{e x + e 3 - e 2 ) - A; 2 (ei + e 2 - e 3 ) 
+ \{ei + e 2 + e 3 ). (100) 



The algebra for the operators (99), (|100|) is found in the work [] I 
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